Manipulation of two spin qubits in a double quantum dot using an electric field 
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We propose purely electric manipulation of spin qubits by means of the spin-orbit interaction (SOI) without 
magnetic field or magnets in a double quantum dot. All the unitary transformations can be constructed by the 
time-dependent Dzyaloshinsky-Moriya interaction between the two spins, which arises from the Rashba SOI 
modulated by electric field. As a few demonstrations, we study both analytically and numerically the three 
operations, i.e., (A) the spin initialization, (B) the two-spin rotation in the opposite directions, and (C) the two- 
spin rotation in the same direction. The effects of the relaxation and the feasibility of this proposal are also 
discussed. 

PACS numbers: 03.67.Lx,73.21.La,71.70.Ej 



I. INTRODUCTION 

Manipulation of spins in semiconductors is a subject of ex- 
tensive studies both theoretically and experimentally. Espe- 
cially, the possible application of the electron spins to the 
quantum computations attracts much attention, and the con- 
trol of a single spin or two spins in quantum dot systems aim- 
ing at the qubit operations with large-scale integration is an 
important issue. Loss and DiVincenzc' proposed the imple- 
mentation of the universal set of quantum gates by using the 
time-dependent exchange interaction and local magnetic field. 
It is known that the XOR and the single-spin operations are 
enough to construct any quantum computations i^ The SWAP 
operation C/swap and its square root U^^i^j^p can be realized 
by the exchange interaction with a certain period of time, and 

1/2 

the quantum XOR gate by the combination of t/gwAP ^^'^ ^^^ 
single-spin rotations induced by magnetic field or a ferromag- 
net. Another proposal is the electron spin resonance transis- 
tors in Si-Ge with the g-factor modulated by the electric field 
serving the possible qubit systemj^ 

Experimentally, the gate voltage can control electron spins 
of a double quantum dot, which contains initialization, ma- 
nipulation, and read-out*l With the help of magnetic field, the 
singlet and one of the triplet form the two-level system, in 
which the SWAP operation and the singlet-triplet spin echo 
have been demonstrated."* 

However, it is desirable to control spins purely electri- 
cally since it is difficult to apply magnetic field confined in 
a nanoscale region. It was proposed to use the decoherence- 
free subspace, in which the exchange interaction alone is uni- 
versal.— Soon later, it was proposed that a single qubit can be 
encoded by three spins, where only two in eight (= 2'^) quan- 
tum states are usedi^ In this proposal, the global magnetic field 
is inevitable in initialization. 

In this paper, we propose fully electric manipulation of 
spins in a double quantum dot, in which magnetic field is not 
necessary at all even in initialization. Furthermore, in our 
method, two-bit operation is realized only by using a dou- 
ble quantum dot in contrast to the previous proposal;^ We 
explicitly show the universal set of quantum gates can be 



FIG. 1 . The schematic (a) top and (b) sectional views of the vertical 
double quantum dot we propose. Shaded and gray regions indicate 
quantum dots and electrodes. Four side gates are attached to each dot, 
which makes it possible to apply electric field in three directions. 



constructed by the exchange and the Dzyaloshinsky-Moriya 
(DM) interactionsii^ 

Most of the electric manipulation methods of spins employ 
the relativistic spin-orbit interaction (SOTjr-^— and it has been 
already demonstrated that the Rashba SOI can be controlled 
by the gate voltage in GaAs systemii^ The Rashba interaction 
is written as 

Hb, — \p- sxE 

with p, s being the momentum and spin of an electron, respec- 
tively, while E is electric field and A is the Rashba coupling 
constant. In a double quantum dot, the Rashba SOI in the re- 
gion between the two dots leads to the spin rotation associated 
with the transfer of the electron, i.e.. 
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C2 + H.C., 



where cj ~ (c. 



^ ct ) 



is the creation operator of the elec- 



tron at ith dot, t is the transfer integral, and e'®''*/-^ is the 
SU(2) matrix corresponding to the spin rotation around the 
axis 9 II 6*12 X E with ei2 the unit vector connecting the sites 
1 and 2. When spin 1/2 is localized at each dot, the transfer in- 
tegral together with the Coulomb charging energy leads to the 
DM interaction [with D \\ 9 in Eq. ([T])] as discussed in semi- 
conductor nanostructures Jiii^ Note that the DM interaction is 
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first-order perturbation in H', it is given by 



FIG. 2. Schematics of the two-qubit operations. (A) Preparation 
of fully polarized triplet state as the initialization from the singlet 
ground state, (B) the two-spin rotation in the opposite directions, and 
(C) the two-spin rotation in the same direction. Black circles indi- 
cate quantum dots, and the black broken oval indicates the singlet 
entanglement due to the exchange interaction J. 



the two-spin interaction, and it appears difficult to manipulate 
only a single spin with use of the DM interaction since the 
single-spin Hamiltonian necessarily breaks the time-reversal 
symmetry. In the ingenious method using the DM interaction 
and the spin anisotropy, the two-spin encoding scheme was 
adopted, i.e., two of four states in every nearest-neighboring 
pair of spins are used. - Also it is shown that the SOI together 
with only one component of magnetic field can manipulate a 
single-spin qubitAiS In contrast, we will show below that it 
is possible to construct the single-spin operations using the 
time-dependent DM interaction without any magnetic field or 
magnetic anisotropy. 



II. PERTURBATION THEORY 

Below we explicitly construct the unitary transformations 
from the exchange and DM interactions, 

H{t) = Ho + H'{t) = J(si-S2-l/4) + 5(i)-si XS2. (1) 

Here the exchange interaction J is assumed to be constant, 
and the DM vector D{t) to have three components, which 
breaks the overall axial symmetry. For the realization of 
this Hamiltonian, we propose the vertical double quantum dot 
schematically shown in Fig. [T] With this configuration, four 
side gates attached to each dot make it possible to produce the 
Rashba SOI and to shift the centers of wave functions in the 
dots separately. The former leads to the DM interaction, and 
the latter is essential to make the z component of the DM inter- 
action.— These two issues are both determined by the spatial 
symmetry, i.e. the original C4 and some mirror symmetries 
with respect the xy plane is lowered by the electric field. We 
show that temporal changes of D{t) enable (A) to initialize 
spins from the singlet ground state to one of the triplet states, 
(B) to rotate two spins in the opposite directions, and (C) to 
rotate two spins in the same direction, which are schematically 
shown in Fig.|2] Combining (B) and (C), we can construct the 
single-spin rotation. 

First of all, we investigate the time-evolution operator U{t) 
within the perturbation theory in the DM interaction since it 
is usually smaller than the exchange coupling J. Up to the 
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e-'""* exp 



dhH{ih 



(2) 



in which H{{t) == e^"'>*H'{t)e~''""* denotes the perturbation 
Hamiltonian in the interaction picture, and we put fi = 1. 
Note that the time-ordered product T in front of the exponen- 
tial operator is absent in the first-order approximation. To ob- 
tain H[{t), we solve the equation of motion. Two spins obey 
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leading to 
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Sll(t) X S2i(t)] 



-Jsii(i) X S2i(i), 
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sii{t) - S2i(t)] = - Jsii(i) X S2i(i). 



2di 
Hence we can explicitly obtain 



H[{t) =D{t) 
1 



ijt 



si X 5*2 COS Jt -\ — (si — 3*2) sin Jt 
2'^/^D-{t)\0){l\+iD'{t)\0){2\ 



_2-i/2^+(i)|0)(3| 



H.c. 



(3) 
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where i:>±(t) = Dy{t)±iD='{t), and \0) = (|t;) - |;t>)/\/2 
is the singlet state, while |1) = jtt), |2) = {\U) + \it))/V2, 
and |3) = \ll) are the triplet states. It is noted here that the 
matrix elements of the DM interaction connect the singlet and 
a linear combination of the triplet states. Considering that the 
Hilbert space of the triplet states has four real degrees of free- 
dom (= 2x3—2 corresponding to the normalization condition 
and the overall phase factor), the three real coefficients D of 
the DM interaction appear to be not enough. This is true if we 
assume that D is independent of time, but once we design the 
time dependence of D{t), we can connect the singlet |0) with 
any linear combination of |1), |2), and |3). 

Let us consider initialization |0) -^ |1). For this purpose 
we take Dini(t) = -D(— sin Jt, cos Jt, 0) and put to = 27r/ J 
in Eq. ^, which makes the first factor e~^^°* in Eq. Q unity, 
i.e., e*'-'*" = 1 leading to 



C/ini = [C/^i^io)]" = exp 



TUT ]~) 

_,^(|0)(1| + |1)(0|) 



\/2J 



When we start from the singlet state, this is a simple two-level 
problem, 

TT ,„, rnrD , , rnrD, , 

C^ini|0)=cos^|0)-zsm-^|l). 

Integer n is determined by the condition mrD/y/2J = tt/2, 
which gives — i|l). This is asymptotically exact in the limit 
of D/J -^ and n -^ 00, and we show in Figs. [3ja) and 
|3jb) the numerical results for finite n by taking all the higher 
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FIG. 3. (Color online) The left panels show time evolution of the spin 
expectation values, in which we choose n = 4 and D/ J — 0.177 for 
(a), n = 4 and D/J = 0.25 for (c), and n = 16 and D/J = 0.199 
for (e). The right ones show n dependence of the fidelity errors de- 
fined in Eq. ([S}. The top two panels (a) and (b) represent initial- 
ization, the middle ones, (c) and (d) represent the two-spin rotation 
around the x axis by ±7r/2, leading to (si/2) |[ ±y> and the bottom 
ones (e) and (f) represent the two-spin rotations around the x axis 
by 7r/2 corresponding to (si/2) || y- The relaxation effect is not 
included. 



-DoppW = (-£> sin Jt, 0,0) in Eq.®. Then we get 



u:^J9) = [u^'Hto)Y 



exp 



' 2 J ^"^1 '^> 



(6) 



which rotates two spins around the x axis by ±6' ~ 
ztmrD/2J. The results of the numerical simulation for fi- 
nite n with 12710/2 J fixed at tt/2 are shown in Figs.[3lc) and 
md). Again the fidelity error is less than 10^^ even for n = 4, 
which decreases as oc l/n'^. 

The two-spin rotation in the same direction is the most non- 
trivial. We define the unitary transformations 



f/"(io)=exp 



■kD 

i — {2{si X S2)"cosv3 - (s" - S2)smip) 



achieved by D"{t) = — Di" cos{Jt + ip) (x" = x, y, z) up 
to the first order in D/J. When we take the "magic angle" 
ip ~ 7r/6, we obtain the following composite operator from 

t/"'s 



c/^: 



M = [u'-uyu-'u-y] 



:CXp 



( 



ttD 



V 2J 



(sT + 4) 



(7) 



in which we apply the Baker-Campbell-Hausdorff formula 



iXAiXB „-iXA„-i\B 



^cxp[-X^[A,B]+0{X^)] 



This is two-spin rotation around the x axis by 6* := 
n(7rD/2J)^. The rotation around the y and z axes can be 
obtained by cyclic permutation. In deriving Eq. (|7]i, it is nec- 
essary that the second-order correction cancels. Up to the 
second-order perturbation, Eq. ^ is modified to 



order terms solving the equation of motion of the density ma- 
trix p{t). The quantum states show oscillatory behavior in the 
time scale corresponding to the exchange interaction J, and 
approach to the predicted states from the first-order perturba- 
tion as n increases, in other words, as D/J decreases. It is 
shown that the fidelity error defined as 



C/(2)(t)=e-»^''*cxp 



1 - ti- 



es) 



is already of the order of 10^'^ for n = 4. Here p is the calcu- 
lated density matrix of the final state, i.e., t = nto = \/2Tr/D, 
while Poo is that of the desired pure state shown in the right- 
hand side of the arrows in Fig. [2] As seen from Fig. Ob), the 
fidelity error is proportional to 1/n^ as expected. This means 
that the conversion between the singlet and triplet can be ma- 
nipulated very effectively by the DM interaction. 

Next we construct the two-spin rotation in the opposite di- 
rections. Let us consider the rotation around the x axis since 
the same applies to that around the y and z axes. Here we take 



dhHiih) 



dti / dt2[H{{h),H{(t2)] 

"^ Jo Jo 
The second-order correction, which is proportional to 



27T pUi 

dui / du2 cos(ui + (p) cos(u2 + f) sin(Mi — U2) 



= -(l-2cos2¥.), 

exactly vanishes at 1^ = 7r/6. As the result, Eq. O holds up 
to the second-order perturbation. Namely, it is asymptotically 
exact in the limit of D/J -^ and n ^^ 00 with n{D/J)^ 
being fixed finite. Figures [3le) andlUf) show the numerical 
results for the finite n cases. The fidelity error in Fig.[3jf) de- 
creases slowly as oc l/n, and therefore the condition is much 
more stringent than that for Figs. |3jb) and[3jb(d). 

By combining Eqs. ^ and (|7]i, we can implement any 
single-spin rotations. In Fig. |4] we demonstrate rotation of 
spin 1 around the x axis by tt/2 with use of [7^(71/4) and 
^samc(^/4)- The former needs the x component of the DM 
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FIG. 4. (Color online) The time evolution of the spin expectation 
values in rotation of spin 1 around the x axis by 7r/2. We set 6 — 
miD/2J = 7r/4 with n = 4 for the operation (B), and then 6 — 
n(7rD/2 J)^ — 7r/4 with n = 4 for the operation (C). Relaxation is 
not included. 



vector, while the latter needs the y and z components. Thus 
three components are required to implement single-spin rota- 
tions around one axis. 

The read-out process can be also designed by the DM in- 
teraction, which connects the singlet and triplet states. This 
is achieved by the transformation from (A'^i, N2) = (1,1) to 
(0, 2) to measure the singlet probability Pq. Here Ni = 0,1 
be the electron number at ith dot. One can always measure 
the probability of singlet Pq by shifting the gate voltage to 
transform (A^i,iV2) = (1, 1) to (0,2)."* By tuning the time- 
dependence of D{t), one can make the situation that the sin- 
glet |0) is coupled to an arbitrary linear combination of the 
triplet states, called |T). The time evolution within this two- 
dimensional Hilbert space spanned by |0) and |T) can ex- 
change the component of |0) and |T) by some rotation of the 
angle tt. From Pq after this rotation, one can measure the 
probability Pt of I T). 



After the Born-Markov approximation and dropping irrele- 
vant terms, we obtain the equation of motion of the density 
matrix as 



dpijt) 
di 



= -rpi(i) 



duti] 



bath 



'Hintlit), [Hintlit - u), pi{t) (g) Pbathi]] 



or 



m^-^[H{t),pit)]-rp{t) 

in the Schrodinger picture. To evaluate Tp, we use 

/•OO 

/ dwtrbath bf{t)b^{t - 7i)pbathe-'"" 
= XI diangjpm / du trbath(a|Q^„e'"'°"* + aia„e"*'^'°"*) 



X (aL^e*'^^'^" (*-")+ a,fl™e-''"^^" 



^j,3m^ 



(t — u)\ —iuju 

^ ')Phs.t\fi 



/•CJO 

=5ij5al3'S^gian^ / du 



X [77.(a;j„„)e*'^'""" + (1 + n{uj^an))e~ 

n 

X [n{Lj.ian)5{Lo - Wjan) + (1 + n(wjQ„))(5(a; + w,a„)] 
=5ij5al3[n{uj)Aia{uj) + (1 + n{-uj))Aia{-uj)], 

in which we define the boson distribution function n{uj) = 
(^g/3w _ 2^-j-i ^jjjj j.jjg inverse temperature /3, and 



III. EFFECT OF RELAXATION 

We choose the integer n to satisfy n{D/J) = const, for 
(A) and (B), and n{D/jy = const, for (C) to set angles. 
There is a trade-off because the DM interaction D/J is as- 
sumed to be small to validate the perturbation theory, while 
we cannot set larger nto than the relaxation time. To substan- 
tiate this consideration by explicit calculation, we employ the 
standard boson-bath model to the relaxation. The spin-boson 
interaction is described as 



in which 



Hint = A(6i • Si +62 • S2), 



bi — 7 ^ gian [O-ian + 0,ian) 



is a fluctuating field from nuclear spins. Here aian is a bosonic 
operator whose motion obeys the Hamiltonian 



is the bath spectral function. After straightforward calcula- 
tions, we get 

(f p)oo =A^[3niAipoo - (1 + ni)Ai{pii + P22 + A33)] 

(f p)ll/33 =A^[(1 + »^l)^lPll/33 - niAipoo] 

+ A^(l + 2no)Ao(pii/33 - P22) 
(f (0)22 =A^[(1 + ni)Aip22 - ni^ipoo] 

+ A2(1 + 2no)Ao(2p22 - Pii - P33) 
(fp)oi_03 =A'[(l/2 + 2?ii)Ai + (1 + 2no)^o]poi-03 
(rp)i2/23 =A2[(1 + ni)Ai + 2(1 + 2no)Ao]pi2/23 
- A^(l + 2no)AoP23/i2 
(fp)i3 -A2[(1 + ni)Ai + 3(1 + 2no)Ao]pi3 

in the singlet-triplet basis. Here uq — n(0) and ni — n{J) 
while ^0 ~ A{0) and Ai = A{J). For simplicity, we ne- 
glect the directional dependence of the spectral function, i.e., 

Aia = A. 

Figure [5] shows the numerical results similar to Fig. [5] but 
with the relaxation. In Figs. |5tb), |3d), and|5lf), the fidelity 
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FIG. 5. (Color online) The relaxation effect is taken into account 
for {2JTi)-^ = (JT2)"^ = 1.0 X 10"'*. The used parameters 
are n = 8 and D/J = 8.84 x 10"^ for (a), n = 11 and D/J = 
9.09 X 10"^ for (c), and n = 18 and D/J = 0.188 for (e), which 
give the minimum fidelity errors, respectively. See the caption of 
Fig.[3]for the notation. 



errors have minimum at a certain n because the contribution 
from relaxation is oc n while the discretization error decreases 
as ex l/n^ or ex 1/n. Especially the two-spin rotation in the 
same direction is greatly affected since it contains four unitary 
transformations and its fidelity error due to the higher-order 
perturbation decreases as 0{n~^). To neglect the effect of 
relaxation, the relaxation time T must be longer than 4nio ~ 
(J/D)^ X (1/J), leading to JT > {J/Df. 



IV. DISCUSSION 

Finally we discuss the realistic setup of our proposal in 
semiconductors. The most serious problem is the relaxation 
of electron spins. The origin of relaxation is mainly the hy- 
perfine interaction with nuclear spins. The effective field 
h^cBnuc is typically SOucV (Ref. 'Ys) and the relative mag- 
nitude J/h'-fcBnuc is a crucial parameter to control the re- 



laxation time. In a GaAs double quantum dot, this value 
is 3-10, and the singlet correlation decays on a time scale 
10ns.— On the other hand, in the case of the singlet-triplet 
relaxation time of the two-electron system in a single quan- 
tum dot, where the splitting is about ImcV and much larger 
than ?i7oi?nuc, the hyperfine interaction is not effective, and 
hence the Ti ^ T2 ^ 0.2msji^ Therefore it is essential to in- 
crease J and to increase Ti, T2. An encouraging theoretical 
analysis gives an estimate of J ^ ImcV in coupled quantum 
dots,^" and also the vertical quantum dots as shown in Fig.[T] 
might enhance J compared with the horizontal dots. 

Another problem is the order of magnitude of the DM in- 
teraction generated by the electric field. According to Nitta 
et al.r^ the Rashba constant can be modified from A = 
0.64 X 10~"eVm at the gate voltage V^ = 1.5V to A = 
0.92 X 10""eVm at V^ = -l.OV. This suggests that it can 
be modified by 1.12meVnmper IV. With the typical distance 
between two dots used in Ref. |23, 2d = 28nm, the SOI is es- 
timated as A = 2(7r/2d) x 1.12 = 0.25mcV at V^ = IV. 
The ratio D/J ^ A/f can reach the order of 0.1. 

Under strong electric field, we may woiTy about the break- 
down phenomenon. However, charge transfer between two 
dots hardly occurs except at the resonance since the energy 
levels of dots are discrete. In addition, the typical time scale 
of charge transfer is of the order of lOOiis according to Ref A 
which is much longer than that of the operations we discuss, 
nto ^ 0.5ns. Therefore we expect that charge transfer can be 
neglected. 

In summary, we have proposed the purely electric manip- 
ulation of qubits in the double quantum dots in terms of the 
Rashba and DM interactions which are modulated by the time- 
dependent voltages. This idea might be useful also for the 
control of macroscopic magnetization in the dilute magnetic 
semiconductor, which is an issue left for future investigations. 
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